Polyhedral colloidal 'rocks': low-dimensional networks 
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We introduce a model system of anisotropic colloidal 'rocks'. Due to their shape, the bonding in- 
troduced via non- absorbing polymers is profoundly different from spherical particles: bonds between 
rocks are rigid against rotation, leading to strong frustration. We develop a geometric model which 
captures the essence of the rocks. Experiments and simulations show that the colloid geometry leads 
to structures of low fractal dimension. This is in stark contrast to gels of spheres, whose rigidity 
results from locally dense regions. At high density the rocks form a quasi one-component glass. 
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Dispersions of mesoscopic colloidal particles are impor- 
tant for several reasons. They model atomic and molecu- 
lar systems, yet single particle level imaging reveals local 
phenomena, such as pinpointing mechanisms of dynamic 
arrest [TJ [2] , which are inaccessible in conventional sys- 
tems. Furthermore, colloidal and nanoparticle systems 
are important materials. Colloidal gels stabilize a range 
of industrial and consumer products from pesticides to 
cosmetics. Underlying both is the tuneability of inter- 
actions between colloids which are theoretically well un- 
derstood and hence enable design of self-assembled struc- 
tures at small lengthscales. Recently, anisotropic interac- 
tions using a number of sticky patches per particle have 
been introduced, opening new routes of self-assembly. 
Reducing the number of patches per particle leads to net- 
works of low fractal dimension [3] . 

Here we consider depletion induced gels of anisotropic 
particles. The addition of polymer mediates an effective 
attraction between the colloids due to excluded volume 
effects - see Fig[TJb). The depth of this attraction is pro- 
portional to the polymer concentration c p and the range 
is set by the polymer size. At sufficient strengths of the 
depletion attraction, spinodal phase separation leads to 
colloid-rich (polymer-poor) 'colloidal liquid' and colloid- 
poor (polymer-rich) 'colloidal gas' phases. Decreasing 
the range of the depletion interaction by reducing the 
size of the polymer relative to the colloid leads to a higher 
density colloidal liquid. If the packing fraction of the col- 
loidal liquid is high enough (0 « 0.58) phase separation is 
arrested. The resulting network of voids combined with 
'arms' of high local colloid density, is termed a gel [4 . 

A range of anisotropic particles has been synthesized 
[5j [6] , and systems with hard interactions have recieved 
considerable attention very recently: rods form the ex- 
pected liquid crystalline phases [7] , anisotropy suppresses 
long-ranged order close to a wall [8] and at high den- 
sity, colloidal dumbells show multiple relaxation path- 
ways [9]. However, although networks of anisotropic par- 
ticles are important, e.g. in shiny paper coatings, sys- 
tematic experimental studies of the effect of anisotropy 
in systems with attractions are limited: attractions be- 
tween colloidal platelets can lead to stacking [10J[TT], gels 
of colloidal rods form bundles [12] and due to their patchy 



interactions, clay platelets form low density colloidal 'liq- 
uids' [13]. 

Here we introduce an anisotropic colloidal system of 
faceted polyhedra, or 'rocks', with tuneable interactions, 
amenable to 3D single-particle level analysis with confo- 
cal microscopy. Unlike gels of spherical particles where 
phase separation is suppressed by slow dynamics due to 
the high local colloid density of the 'arms', we show that 
the polyhedral nature of the rocks leads to bonds which 
do not rotate and thus rigid structures and networks of 
low fractal dimension are formed. It appears that for this 
geometry, a modest depletion attraction can drive the 
system into a regime similar to diffusion-limited cluster 
aggregation (DLCA) where the bonding is irreversible. 
Inspired by the remarkable effect of particle geometry, 
we develop a model which captures the essential geomet- 
rical properties of the colloidal rocks in the form of do- 
decahedral clusters of 20 spheres, as shown in Fig. []Je). 
We further explore the effect of particle geometry at high 
density and find that the colloidal rocks are a quasi one- 
component glass- former. 



I. EXPERIMENTAL 

Faceted polytetrahedral rock-shaped particles as 
shown in Fig. [TJa) were produced by modifying a synthe- 
sis to produce poly-methyl-methacrylate (PMMA) col- 
loids stablised with poly- hydroxy stearic acid [H] . De- 
tails of the synthesis will be presented elsewhere. The 
particle size is expressed as a r = 3.5 /im, the longest 
dimension of each rock in images such as Fig. [TJa). It 
is hard to determine the polydispersity for these parti- 
cles, which have both shape and some size polydispersity. 
However, inspection of SEM images suggests that the size 
polydispersity is ~ 5%. The colloidal rocks were labelled 
with rhodamine iso-thiocynate fluorescent dye. To in- 
vestigate the change induced by anisotropy we also used 
spherical PMMA particles, of diameter a s = 2.4 /im with 
4% polydispersity determined by static light scattering. 

The colloids were dispersed in a density- and refractive 
index matching mixture of cis-decalin and cyclohexyl bro- 
mide (CHB), to which was added 4 mMol of tetrabutyl 



Figure 1: (color online) (a) SEM image of colloidal rocks, (b) Schematic of the mechanism of rigid bonding induced by multiple 
overlap zones in the rocks contrasted with spheres which are free to rotate. Regions from which the polymer coil center of 
mass is excluded are shown in pink, (c) Confocal image of a gel of rocks c p /cp el — 1.13. (d) Confocal image of a gel of spheres 
c p/ c p el — 1-44. c p is polymer concentration and c g v el is at the onset of gelation, see Fig. pi (e) Dodecahedral assembly of 20 
spheres used in simulations. Scale bars in a,c,d are 10 /im. 



ammonium bromide salt to screen electrostatic interac- 
tions [2 . Polystyrene polymer was added to induce deple- 
tion attractions. The polymer used was M w = 2.1 x 10 T 
g/mo in the case of the rocks and M w = 3.1 x 10 T 
g/mol for the spheres resulting in polymer radii of gy- 
ration of R g = 180 nm and R g = 220 nm [15] and a 
polymer-colloid size ratio 2R g /<j of 0.11 and 0.18, re- 
spectively. With the exception of the modest change in 
interaction range, the experimental conditions were iden- 
tical for both rocks (r) and spheres (s). The colloid vol- 
ume fraction <p r = <j) s = 0.05 was determined by weighing 
out the samples. We track the colloid coordinates in 3D 
using confocal microscopy images. The rocks are suffi- 
ciently large with an intensity maximum at the centre 
of the particle that algorithms developed for spheres [16] 
identify the rock centres effectively with an error < 100 
nm. Particles were defined as bonded if their centres lay 
within the interaction range a + 2R gi i.e. 3.9 /im and 
2.8 /im for the rocks and spheres respectively. Moderate 
changes in this bond length had no significant impact on 
our results. 

II. RESULTS AND DISCUSSION 

Geometry introduces two key differences between the 
rocks and spheres relevant to depletion attractions as in- 
dicated in Fig. [IJb). Firstly, for rocks the overlap of 
excluded volume, which drives the depletion interaction, 
is much reduced compared to spheres. This is due to the 
concave shape of the rocks which means that much of 
the volume from which the polymers are excluded can- 
not overlap due to the approach of two rocks. This is 
shown by the pink regions in Fig. [ljb) which denote ex- 
cluded volume, while the black is the colloidal particle. 
Blue denotes regions where overlap of excluded volume 
is prevented due to the concave faces. Thus we expect 
that more polymer is required in the case of the rocks, 
for a given degree of attraction. Secondly, the faces of 
the rocks lead to multiple overlap zones as shown by the 
overlapping pink regions in Fig. [jjb). These impose an 



energetic penalty to rotation, which is absent in the case 
of spheres. Both points were confirmed by calculations 
with model rocks and ideal polymer. 

Confocal microscopy images of gels of rocks and 
spheres are shown in Fig. [TJc) and (d). The difference 
is striking: rocks form one-particle wide chains while 
spheres form densely packed structures. This provides 
direct evidence for a different bonding scenario in the 
two species, as indicated in Fig. [ljb) . We identify gela- 
tion with a percolating network of long-lived bonds. Gel 
dynamics are shown in supplementary movie 1. In fi- 
nite sized microscope images, it can be hard to directly 
measure percolation. We therefore determine percolation 
through the ratio L c /L as shown in Fig|2|a), where L c is 
the length of the largest connected cluster and L is the 
image size. Gelation is thus identified with a rapid rise 
in the ratio L c /L. The polymer mass fraction for gela- 
tion c g v el in the case of spheres is 4.1 ± 0.5 x 10 -4 and 
2.0 ± 0.3 x 10 -3 for the rocks, which corresponds to a 
volume fraction (j)^ el of 0.42 ± 0.05 for the spheres and 
2.3 ± 0.2 for the rocks. Here the polymer volume frac- 
tion is defined <p p = Aitp p R^ g /A where p p is polymer num- 
ber density. We attribute this fivefold change in poly- 
mer volume fraction required for gelation to the much 
reduced overlap volume in the case of the rocks. Effects 
due to different polymer-colloid size ratios and polymer 
non-ideality were estimated to be small [T7] . 

The difference between the rocks and spheres in the 
structure at the particle level is also striking, as shown 
by the radial distribution functions in Fig. |3ja) and 
(b). Data for spheres has been compared against com- 
puter simulation for hard spheres in [15], showing excel- 
lent agreement. The rocks show similar behaviour, in- 
dicating that the colloids are stable against aggregation. 
At higher polymer concentration, although both species 
show a strong increase in the first peak upon gelation, in 
the case of spheres, the onset of gelation (c p /c 9 p el > 1) 
coincides with g(r) > 1 for r < 5<r s , indicating clus- 
tering/dense regions in agreement with Fig. |TJd) . This 
behaviour is entirely absent in the case of the rocks [Fig. 
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Figure 2: (color online) Percolation and fractal dimension, (a) Size of largest connected region L c scaled by image size L. 
Dashed line is a guide to the eye. (b) Radius of gyration of clusters as a function of the number of particles in the cluster 
plotted to yield the fractal dimension df (experimental data from rocks), df in experiments (c) and in simulations (d) for 
isolated clusters with <j) r — 4>s — 0.0125. In (a,c,d) the rocks are denoted by pink squares □ and the spheres by circles •. 



[TJc)]. In other words, for spheres, phase separation has 
already begun before gelation sets in: gelation is driven 
by the arrest due to the development of the dense phase 
[4], once it has formed on a lengthscale of several parti- 
cles. The rocks show no signs of locally dense regions: 
gelation is driven soley by rigid bonds. Similar conclu- 
sions may be reached by considering the number of neigh- 
bours of a given particle [Fig. [3jc,d)]. In the case of 
spheres, the onset of phase separation (gelation) leads to 
a considerable change in the number of neighbours, in- 
dicating locally dense regions. Rocks on the other hand 
show no such indication of phase separation. In other 
words, rocks have a reduced valency induced by their ge- 
ometry with some similarity to patchy particles [3]. 

This bonding and local structure has significant im- 
plications for the network formation. To quantify the 
structure, we use a local 'fractal dimension', df, for finite 
assemblies of particles, which is defined as R c g oc N x ^ df 
where R c g is the radius of gyration of a connected region 
of N particles, as shown in Fig. |5|b). At the low pack- 
ing fractions we consider, the fluid phase forms clusters, 
whose df is found in the same way. We plot the result- 
ing df in Fig. [2|c). Spheres and rocks show an increase 
in df upon increasing the polymer concentration, up to 
gelation, after which non-ergodicity leads to a failure to 
relax locally and a decrease in df, as found previously 
[18] . However, the fractal dimension of rock clusters and 
gels is significantly lower than that of spheres. This is 
consistent with the idea that the sphere-based clusters 
have already begun to condense (leading to an increase 
in df), while the rock clusters show no signs of such 
condensation. We note that for the rock gels, the mea- 
sured fractal dimension is even less than that expected 
for DLCA where df = 1.8. This is surprising: DLCA as- 
sumes unbreakable bonds, which is entirely different from 
the weak depletion-induced gelation expected in colloid- 
polymer mixtures. 

We have argued that the only relevant difference be- 
tween the rocks and spheres lies in their shapes and now 
enquire as to the geometrical origins of the rock be- 
haviour. Unlike other approaches to produce networks 
with low fractal dimension which rely on a few sticky 
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Figure 3: (color online) Radial distribution functions g(r) for 
rocks (a) and spheres (b), respectively. Dashed line in (b) 
denotes the approximate lengthscale of incipient phase sepa- 
ration for the spheres. Data are shifted for clarity. Distribu- 
tions of the number of neighbours P(N n ) for rocks (c) and 
spheres (d), respectively. Labels denote c p /cp el . 



patches on each particle to reduce the number of bonds 
[3], here we consider only the particle geometry with no 
restriction on the number of neighbours a particle can 
have, except the steric limitation which also applies to 
particle with a spherically symmetric interaction. We 
introduce a model based on dodecahedral assemblies of 
20 overlapping spheres whose centres lie on a sphere, as 
shown in Fig. [IJe). This model was choosen as a sim- 
ple representation of the rock shape, as it is reasonably 
resistant to ordered packing, and the faces should mimic 
the rigid bonding of the rocks. These (small) spheres 
have a diameter of 2/7a c , where a c is the diameter of the 
circumscribed sphere. The volume of a set of overlaping 
spheres can be calculated exactly by e.g. the Connolly 
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algorithm [19] and is approximately 0.416 o\ for these 
model rocks. 

For short ranged attractions, the detailed form of the 
interaction potential is irrelevant for many properties 
[20] , so we treat the attractions that result from the ad- 
dition of polymer with a square well of width 0.042<j c 
and depth e. We also carry out simulations of spheres 
of diameter a s interacting via the square well potential 
of width 0.042a s . Using the extended law of correspond- 
ing states, [20], this maps to a polymer-colloid size ratio 
of 1/7, comparable to the experiments. Here we focus 
on the geometric properties, rather than the dynamics: 
our purpose is to demonstrate the qualitatively different 
nature of aggregation of spheres and rocks. We there- 
fore use Monte Carlo (MC) simulations in the canoni- 
cal ensemble with 512 particles. Simulations with more 
(2048) particles showed minimal differences, as for the 
state points we consider. We use a short (0.04<j) and 
long (cr) translational step for both spheres and rocks. 
The former enables local restructuring of bonded parti- 
cles, the latter accelerates aggregation. For the rocks, 
small (0.01 rad) and large (0.2 rad) rotational moves are 
also used. Each simulation is 'equilibrated' for 10 4 MC 
sweeps prior to sampling. Each sweep involves attempted 
moves of all types. Simulations are repeated at least six 
times. 

In Fig. 2^d) we plot fractal dimensions for simula- 
tions of rocks and spheres for packing fraction c/> r = <j) s = 
0.0125 (the experiments are carried out at c6 r = (j) s = 
0.05) [21 . The spheres exhibit similar behaviour as in 
the experiments, with df ~ 2 when the system falls out 
of equilibrium on the simulation timescale. This we take 
as the strength of attraction associated with spinodal de- 
composition [4], which corresponds to fie* « 2.53, noting 
that for such short-ranged interactions, the spinodal line 
depends rather weakly on c/>. For the rocks, we identify 
e* & 1.75 fc#T with the onset of a regime where the 
potential energy decreases continuously. Both rock and 
sphere df are lower than the experiments at low attrac- 
tion, which may be related to the change in <j>. 

The simulations clearly show that the rock clusters 
have a fractal dimension which is less than that of the 
spheres [Figs. |5Jc,d)]. In fact df ~ 1.5, even less than 
in the experiments. Although the simulation results are 
presented for a lower density the agreement between ex- 
periment and simulation adds further weight to our hy- 
pothesis that (i) the rocks have an intrinsically differ- 
ent bonding mechanism compared to spheres driven by 
geometry and (ii) that their behaviour is qualitatively 
captured by the dodecahedral clusters. Since the rocks 
can bond to as many neighbours as they have faces, they 
they can form dense phases and therefore, at equilibrium 
dense and dilute phases coexist so the rock gels should 
be metastable. 

The fractal dimension of rock clusters, both in experi- 
ments with percolating networks and simulations of iso- 
lated clusters is low, ~ 1.7 and ~ 1.5 respectively, lower 
even than DLCA. Cates et.al. [22], suggested that sim- 
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Figure 4: (color online) (a) Mapping g(r) of rocks at <j) r — 0.48 
(pink squares) and g(r) of spheres at (j) s = 0.54 (black cir- 
cles) from experiments. cr e ff denotes a r and a s for rocks and 
spheres respectively. We estimate the mapping by compar- 
ing the higher-order peaks. At the time of measurement the 
metastable sphere fluid showed no signs of crystallisation, (b) 
Confocal image of rock glassy state. Bar=10 /im. 



ilar behaviour might be found in spheres, if they were 
suddenly - and deeply - quenched to the regime where 
DLCA dominates. The rock geometry may be inter- 
preted as exhibiting similar behaviour with much weaker 
attractions. 

Finally, we consider the high density behaviour. In 
experiments without polymer, rocks did not crystallise 
at a packing fraction c6 r = 0.46 on a timescale of one 
week (Fig. [4|. This is reasonable, as it has recently been 
shown that only a limited degree of asphericity is required 
to suppress crystallisation [23 . Supplementary movie 2 
shows that rocks at c6 r = 0.46 exhibit glassy dynamics, 
with displacement hindered by their neighbours. We find 
that the radial distribution function of dense amorphous 
packings of rocks at <\> r — 0.46 has peak locations and 
a decay of the oscillations similar to the g(r) of spheres 
at <j) s = 0.54 [Fig. j4^a)] . This indicates that the longer- 
range structure is similar, suggesting an effective packing 
fraction for the rocks higher than c6 r = 0.46, likely due 
to the concave geometry. The failure of the rocks to 
crystallise indicates that they are a 'quasi one-component 
glass former'. 



III. CONCLUSIONS 

We have developed an experimental model system 
of anisotropic colloidal rocks whose interactions can be 
tuned and which can be visualized directly in 3D. Due 
to their shape, the bonding introduced via the deple- 
tion attraction is profoundly different to spheres: the 
bonds formed between rocks are rigid against particle ro- 
tation. The frustration induced by this change in geome- 
try profoundly influences the kinetic pathway: rocks form 
open networks of low fractal dimension. We presume 
that the rock gels are, like gels of spheres, ultimately 
metastable to phase separation, but expect that the rock 
gels should be very long-lived. We demonstrated that ge- 
ometry is the dominant driver for this behaviour by intro- 
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during a model of dodecahedral clusters which captures 
the essence of the colloidal rocks. At higher density the 
experimental system forms a quasi-one-component glass. 

Our work opens up a single particle level approach to 
tackle important problems such as the role of friction 
and surface roughness in colloidal and nanoparticle pro- 
cessing and jamming. We hope to stimulate further ex- 
periments, for example the use of external fields such as 
shear, along with theoretical and simulation work which 
may fully determine the phase diagram, and the role of 
particle geometry in kinetic trapping and its interplay 



with hydrodynamic interactions. Self assembly of nan- 
odevices requires avoidance of kinetic trapping [24 ; here 
we have shown the potential of direct visualisation as 
a means to develop control of self-assembly which may 
yield insight into nano-assembly. 
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